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Abstract 

Using a self-similar approach a general nonsteady theory is elaborated for the case of the diffusion 
growth of a gas bubble in a supersaturated solution of gas in liquid. Due to the fact that the solution 
and the bubble in it are physically isolated, the self-similar approach accounts for the balance of 
the number of gas molecules in the solution and in the bubble that expells incompressible liquid 
solvent while growing. The rate of growth of the bubble radius in its dependence from gas solubility 
and solution supersaturation is obtained. There is a nonsteady effect of rapid increase of the 
rate of bubble growth simultaneous with the growth of the product of gas solubility and solution 
supersaturation. This product is supplied with a limitation from above, which also stipulates 
isothermal conditions of bubble growth. The smallness of gas solubility is not presupposed. 
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Introduction 



The basis of the study of the diffusion growth of a gas bubble in a supersaturated solution 
of gas in liquid is usually l[ [2], [], Q, ^, [f| founded on a steady solution of the diffusion equation 
of gas molecules in liquid. It is the purpose of this paper to elaborate a general nonsteady 
theory for the diffusion growth of a gas bubble in a supersaturated solution of gas in liquid, 
exploiting a self-similar approach elaborated in [?J by the example of the diffusion growth 
of a liquid droplet in a supersaturated vapor-gas mixture. This theory does not exploit 
the supposition that the concentration profile of gas dissolved around the growing bubble is 
steady. 

The object under consideration is a liquid solution of gas which includes a gas bubble 
of the same single- component substance that is dissolved in the liquid. We shall consider 
the volume of the solution to be so large that heterogeneous nucleation of bubbles at its 
boundaries does not influence the volume of the solution where the bubble is located. The 
bubble is rather large, which makes it possible to neglect the Laplace forces. The solution is 
presupposed to be diluted, while the gas in the bubble - perfect. Dissociation and chemical 
transformations of the dissolved molecules are also neglected. Due to high heat conductivity 
of the liquid solution, we shall consider all the system of the solution and of the bubble as 
having the same temperature. Temperature and pressure in the solution are given. The 
solubility of gas in liquid is traditionally seen as a dimensionless value defined (at given 
temperature and pressure) as the ratio between the volume of gas and the volume of liquid 
dissolved in this gas. 

We shall study a realistic situation when the diffusion of molecules of the dissolved gas in 
the presence of a bubble occurs simultaneously with the motion of the liquid solvent which 
is caused by the movement of the bubble surface and due to the solvent incompressibility. 
This problem has not been studied in . Taking into account the fact that 

the solution and the bubble in it are physically isolated, we shall consider the balance of the 
number of gas molecules in the solution and in the bubble that expells the incompressible 
liquid solvent while growing. Using the self-similar theory as a basis, we shall derive an 
equation for the rate of growth of the bubble radius in its dependence on gas solubility 
and solution supersaturation. We shall describe a significantly nonsteady effect of a rapid 
increase of the rate of bubble growth at the increase of the product of gas solubility and 
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solution supersaturation. This product will be provided with the limitation from above, 
which stipulates the isothermal conditions of bubble growth. The smallness of gas solubility 
will not be presupposed. 

I. GAS BUBBLE GROWING IN LIQUID SOLUTION OF GAS: GENERAL IDEAS 

The state of the solution is stipulated by temperature T, pressure II and the initial 
concentration (the number density of molecules) no of the dissolved gas. The meaning of 
value n will be revealed further in the initial condition in Eq. (fT2|) and in the boundary 
condition in Eq. ( lT5i) . will stand for the concentration of the dissolved gas in the 
saturated solution, which at given temperature T and pressure II is characterized by chemical 
and mechanical equilibrium with the clean gas above the plane contact surface. We shall 
consider the solution to be diluted. The solution contains a gas bubble (the same gas that 
is dissolved in the solution). Due to high heat conductivity of the liquid solvent, the bubble 
of gas will have the same temperature T as the solution does. The condition that the latter 
is valid and that the temperature in the bubble is homothermal will be justified in section 
IVl There we shall also justify the mechanical equilibrium of the bubble and the solution. 
The radius of the bubble will be designated as R. 

We shall suppose the radius R to be that large that the following strong inequalities are 
observed: 

#>2(r/n, (i) 

R^>AsD/a c v T . (2) 

Here a is the surface tension of the liquid solvent (in the case of a diluted solution), D is the 
diffusion coefficient of gas molecules in the liquid solvent, a c is the condensation coefficient 
of a gas molecule at the virtual transition of gas molecules from the bubble and into the 
solution vt is a mean heat velocity of a gas molecule inside the bubble, s is the gas 
solubility that is defined as a dimensionless value via 

s = kTrioo/Tl, (3) 

where k is the Boltzmann constant. 

When Eq. ([ID is observed, the influence of the Laplace forces on the bubble is rather small 
and the pressure of gas inside the bubble at its mechanical equilibrium with the solution 
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is equal to pressure II of the solution. At the same time the gas inside the bubble can be 
considered perfect; therefore for gas concentration in the bubble, that is designated as n g , 
we have 

n g = II/fcT, (4) 

so n g is unambiguously defined by temperature T and pressure II of the solution. Using Eq. 
(j3J), we can also rewrite Eq. (j4j) as 

n g = n 00 js. (5) 

Due to Eq. dSJ) value s defined in Eq. (J3J) corresponds to the usual understanding of solubility 
that is, a ratio (at given temperature and pressure) of the volume of gas to the volume of 
liquid that has dissolved this gas. Solubility s and value depend not only on temperature 
T and pressure II, but also on the type of solvent and on the type of dissolved gas. 

We shall neglect the fugacity of the liquid solvent and do so for the following reason: 
let Pp be the pressure of the saturated vapor of the solvent. If we have a look at the 
thermodynamic diagram, at temperature Tand pressure II of the solution the solvent will 
be found deep within the area of the steadiness of the liquid phase. For this reason Pp<gjl, 
which, using concentration np of the saturated vapor of the liquid solvent, can be rewritten 
as rig-Cll/fcT. From this and from Eq. (T4j) it follows that n^-Cn^. Thus the bubble consists 
mainly of gas; that is why the fugacity of the solvent can in fact be neglected. 

As it was shown in [3], the strong inequality in Eq. (T5]) is the condition for the existence 
of the diffusion regime of bubble growth. Leaving a certain reserve, let us assume that 

i2b~20nua(^,— }. (6) 
L II a c v T J 

Eq. ([6]) evaluates value Ro of the bubble radius, which is initial for the existence of conditions 
in both Eq. (DO) and Eq. (j2J). Jointly with those, Eq. (T4|) and the diffusion regime of the 
bubble growth are also well observed. The largest of the two values given in braces in Eq. 
(JHJ) is, as a rule, the value of 2er/IL 

Let us introduce characteristic time to, at which the following initial condition on the 
radius of the bubble is observed: 

R(t)\ t=t0 = R . (7) 

Time t is calculated from the time of nucleation of the bubble that is growing irreversibly. 
From Eq. (JTj) and the meaning of value Ro it can be seen that after time to after the 
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nucleation of the bubble Eq. (j4]) and the diffusion regime of bubble growth will be observed 
simultaneously. The larger Rq, that is, the more precise is the observation of Eqs. ([1]) and 
(T2J), the larger t . 

The bubble perturbs the surrounding solution. Let us designate the concentration of the 
dissolved gas at distance r from the centre of the bubble at time t as n(r,t). If the growth 
regime of the bubble is a diffusion one, the evolution equation is observed 

= DAn (r, t) - div [n(r, t)v(f, *)] (m ). (8) 

The first member in Eq. (jSJ) accounts for the diffusion of gas molecules, while the second 
member stipulates it that this diffusion occurs against the background of the movement of 
liquid solvent caused by the movement of the surface of the bubble and by the incompress- 
ibility of the solvent. If we study the velocity of solvent movement designated in Eq. (jSJ) as 
v(r, t) in the spherical system of co-ordinates with the reference point in the centre of the 
bubble, we shall see that it has only radial component v r (r, i), for which, as it is obvious, 
the following is valid: 

, , R 2 (t)dR(t) 
Using Eq.Q and some well-known formulas of vector operations 



let us rewrite the evolution equation from Eq. (jSJ) as follows 



dn(r,t) D d 
dt r 2 dr 



2 dn(r,t) 



dr 



R 2 (t) dR(t) dn(r,t) 



(m ). (ii) 



r 2 dt dr 

At all r > we have the following initial condition on the homogeneity of the concentra- 
tion profile of the dissolved gas at time t = of the bubble nucleation: 

n(r,t)\ t=0 = n (r > 0). (12) 

This condition brings about value n as the initial concentration of the dissolved gas. In a 
supersaturated solution it is observed that n > n^. 

During the diffusion regime of the bubble growth we have an equilibrium boundary con- 
dition on the surface of the bubble: 



n{r,t)\ r=R(t) = n 00 (t^t ). 
5 



(13) 



We shall neglect the dependence of the right side of Eq. (fl3l) on the radius of the bubble 
when this radius is significant and conforms with the limitation from Eq. ([1]). 

Let us impose a boundary condition of the absence of an external source of the dissolved 
gas: 

r 2 = o. (14) 



dr 

This condition means that the solution and the bubble in it are physically isolated. 
In accordance with Eqs. (|T2|) and (TH|) it is also true that 

n(r,t)\ r=00 = n . (15) 

Due to physical isolation of the solution and the bubble in it, value no also has the meaning 
of concentration of the dissolved gas at an infinite distance from the bubble; and this value 
is constant in time. 

Each element of the surface of the bubble of radius R is moving with the velocity of the 
movement of the incompressible liquid solvent. Thus, 



dR(t) _ dn{r, t) 

Tin T. — i-J 



dt dr 



(t>t ). (16) 

r=R(t) 



Eq. (1161) can be interpreted simply. It corresponds to the diffusion flux of the numbers of 
molecules through the surface of radius R (there is no convection flux due to mobility of this 
surface). The diffusion flux is spent on the change in the number of molecules of gas in the 
bubble (which at t^to , in concordance with the initial condition from Eq. (17!) . has constant 
concentration n g of pure gas). As the sign of the derivative by r in Eq. (1TB1) obviously 
corresponds with the sign of value n — n^, then, in the case of supersaturated solution, Eq. 
(1T6!) leads to 

dR/dt > (t^t„) (17) 

- i. e. leads to the growth of the bubble with time. 

In order to specify the state of the solution one can use value ( instead of value n . Value 
C can be defined as 

C = (no - noo)/noo. (18) 

In a supersaturated solution, where ( > 0, value ( will be referred to as the supersaturation 
of a solution. 
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II. SELF-SIMILARITY OF THE THEORY 



[71 



Following [7j, let us introduce a self-similar variable p having assumed that 

» = W) <19) 

We shall look for the solution n(r,t) of the evolution equation (ITTj) in the form of n(p) of 
one variable p: 

n(r,t)=n(p). (20) 

Here and further, for the sake of simplicity, we shall not mention that p is changing within 
the limits of p^l. 

In accordance with Eqs. f|T9|) and fTSOj) we have 

dn(r,t) r dR(t) dn(p) 



dt R 2 (t) dt dp 

dn(r,t) 1 dn(p) 

dr ~ R(t) dp 

From Eqs. ( fl6l) and ( |22i) we have 

dR(t) D dn(p) 



(21) 
(22) 



dt R(t)n g dp 

Let us rewrite Eq. ff23l) as 



. (m ). (23) 

p=l 



dR(t) Db 
or, equivalently, as 

dR 2 /dt = 2Db (t^t ), (25) 



where we introduce an important dimensionless parameter b using 

1 dn(p) 



n g dp 



(26) 



P =i 

In accordance with Eq. fT25|) . the derivative dR 2 /dt does not depend on time: along the axis 
of variable R 2 the bubble is "moving" with the velocity that does not depend on time and 
on the size of the bubble. 

Relative measure of inaccuracy of the self-similar velocity of the movement of the bubble 
2Db in Eq. (1251) is equal to the relative measure of inaccuracy of Eq. (Tj|]), which, in its turn, 
does not exceed value 2a /RqII, which can be easily obtained from the estimate from Eq. 

(ED. 



It is easy to make sure that, if we take Eqs. (fT9l) - (1221) and Eq. (}!MI) into consideration, 
the evolution Eq. (fTTl) turns into an ordinary differential equation for function n(p) of the 
following type: 



dp 2 



2 J 1 
- + b[p- — 

P V P 



d 4^ = 0. (27) 



Integrating Eq. ( 1271) . obtaining a constant from the first integration using Eq. ( l26l) . and a 
constant from the second integration using the boundary condition from Eq. ({TBI , formulated 
by virtue of Eqs. f[T$j) and (T2"0"j) as n(p)| p=1 = n^, we have 

n(p) = noo + n g be^ 2 f %^ V 2 "^ . (28) 

J i x 

The fact that Eqs. f[2~7"l) and (12^1) refer only to times t^t is not stipulated for the sake of 
simplicity. As it is obvious from Eqs. ffl9l) and (1201) . a self-similar formula (1281) complies 
with the boundary condition from Eq. (fl^l) . that is the absence of an external source of the 
dissolved gas. 

In order to find parameter b using Eq. (1281) and the boundary condition from Eq. (fT5|) 
formulated by virtue of Eqs. ( TL9l) and (1201) as n(p)| p=00 = n , we come up with the following 
transcendental equation: 

dx 



a = 6e 36/2 / ^ e -^/ 2 -V^ (29) 
Ji x 

where an important dimensionless parameter a is introduced using 

a = (30) 

n g 

In the case of a supersaturated solution, from Eq. fl3"Ul) it follows that 

a > 0. (31) 

In accordance with Eqs. ([5]) and (fl8l) . we can rewrite Eq. (|30l) as: 

a = s(. (32) 

In order that the integral in Eq. (1291) should converge and, simultaneously, that the right 
part of Eq. (1291) should be positive in concordance with Eq. (13T1) . it is required that 

b > 0. (33) 

Eqs. (pHl) and (1331) have the previous inequality ([17]) as a consequence. 




FIG. 1: Dependence of the solution b of Eq. (|29p on parameter a. 

Eq. (1291) for b relative to a has an unambiguous solution within the area stipulated by 
Eq. fEJTj) . The results of numerical solution of Eq. (1291 are represented by the curve in Fig. 
1. This curve does not depend on the type of the solvent and on the type of the dissolved 
gas. All initial parameters of the problem are represented on the curve solely by means of 
parameter a introduced in Eq. fl30|) . 



9 



Integrating Eq. (1251) with the initial condition (JTj), we have 



R 2 (t) = 2Db- (t-t ) + Rl (t^t ). 



(34) 



When Rq is significant, which is required by the estimate from Eq. ([6]), the relative measure 
of inaccuracy 2a /RU of the self-similar rate of bubble growth 2Db in Eq. fl25|) does not 
exceed 10% at times t > to/4 (and does not exceed 5% at times t > t ). During the longest 
stage of the whole time interval < t < to, that is to/ 4 < t < t , the bubble, consequently, is 
growing practically with a self-similar rate 2Db. For this reason we can use an approximate 
formula 



where it is stipulated that the bubble is growing with a self-similar rate 2Db throughout the 
whole interval of < t < t . From Eqs. ( J34l) and ( l35l) we have 



which makes it possible to findi?(t) in Eq. (JT9l) at all times t^to- The larger Ro, the more 
accurate Eqs. f[3"5l) and fl3"ET) are. 

Let us now explicitly formulate the sequence in which the values that are found during 
an experiment are introduced within the theory. The state of the solution is given by its 
temperature T, pressure II and the initial concentration n Q of the dissolved gas (introduced 
into the theory by means of the initial condition in Eq. (1.12)). Using the given T and II 
one can obtain from an experiment (at a stipulated type of solvent and type of a dissolved 
gas) value n^, and, knowing T, II and n^, using definition from Eq. ([3]), we introduce the 
theory with gas solubility s. Then, using either of Eqs. (BJ or ([51), we obtain value n g (that 
is unambiguously defined by temperature T and pressure II). Finally, knowing n , and 
n g , by means of definition from Eq. (!30|) we introduce value a. 

It should be mentioned that in order to stipulate value n in the experiment it is conve- 
nient to proceed from saturated solution at pressure that is larger than pressure II that is 
of interest to us. Let us identify value n,^ obtained during the experiment with value n , 
and then instantly reduce the pressure in the solution to II. Due to low compressibility of 
the liquid solvent, the solution will become just what we are interested in - supersaturated 
solution with the given concentration n of the dissolved gas. 



to~Rl/2Db 



(35) 



R 2 {t)~2Dbt (t^t ), 



(36) 
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III. DEPENDENCE OF THE BUBBLE GROWTH RATE ON GAS SOLUBILITY 
AND SOLUTION SUPERS ATURATION 



In accordance with section [III self-similar Eq. (I29p defines parameter b as an unambiguous 
function b = b(a) of parameter a throughout all the area stipulated by Eq. f[3"Tj) . Function 
b = b(a) does not depend on the type of the solvent and on the type of the dissolved gas. 
Having multiplied function b = b(a) by 2D, using self-similar Eq. ( 125]) . we have the following 
equality 



for the constant in time derivative dR 2 /dt, which characterizes the rate of bubble growth 
(the rate of the growth of its surface). Function b = b(a), which was obtained via numerical 
solution of the self-similar Eq. (129]) . is just what is presented in Fig. 1. It can be seen that 
with the increase of parameter a in the area stipulated by Eq. (I3TI) this function is growing 
faster than parameter a is. 

The physical meaning of parameter a is revealed in Eq. (1321) . Exploiting this in Eq. (l3Tj) . 
we have 



Eq. (1381) shows that derivative dR 2 /dt is proportional to diffusion coefficient D of the 
molecules of gas in the solution and depends on the solubility of the gas s and solution 
supersaturation ( only via the product of s(. Curve 1 in Fig. 2 represents data concerning 
the dependence of dimensionless value D~~ x dR 2 jdt on s( obtained by means of Eq. (I3"8j) 
and the numerical solution of the self-similar Eq. (T2T?j) . It is clear that with the increase 
of s( from s( = 0,01 to s( = 20, value D~ x dR 2 jdt is growing at an increasing rate from 
D- 1 dR 2 /dt = 2 • 10~ 2 to D~ 1 dR 2 /dt = 1,5 • 10 3 . The revealed effect of the significant 
increase of the rate of bubble growth simultaneous with the increase in the product of gas 
solubility s and solution supersaturation ( is nonsteady by its virtue; and this will be shown 
in section HVl 

When a>10, the solution of the self-similar Eq. (1291) can be found analytically. Let us 
rewrite this equation as 



dR 2 /dt = 2Db{a) (t^t ) 



(37) 



dR 2 /dt = 2Db(s() (t^t ). 



(38) 




(39) 
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FIG. 2: Dependence of dimensionless value D 1 dR 2 jdt on the product s(. 

Expanding the exponent from Eq. (|39|) into a Taylor series in powers of x — 1 in the 
vicinity of x = 1, noting that at b^>l the exponent "cuts" the integral in Eq. fl39|) already 
before 1/x 2 manages to practically deviate from unit, using saddle-points technique, we have 

f°° dx ( bx 2 b 3b\ /vr\i/2 
Substituting Eq. fj40l) in Eq. fl39|) . we come up with an analytical solution of the self-similar 



Eq. (J29D: 



b(a) = -a 2 (a>10). (41) 

7T 



Eq. ( I4TI) at the limit of its applicability a~10 concords well with the graph given in Fig. 1. 
Thus analytical Eq. (HIT) conjointly with Fig. 1 provides comprehensive data on function 
b = 6(a) throughout all the area stipulated by Eq. ( 13TI) . Curve 2 in Fig. 2 represents data 
on the dependence of dimensionless value D~ l dR 2 /dt on s( obtained by means of Eq. (138]) 
and the analytical Eq. (l4Tj) . 

Exploiting Eq. (j41]) in Eq. (1381) . we have 

^ = — «) 2 «>10, t^o). (42) 
at 7r 

Analytical formula fj42|) shows that the derivative dR 2 /dt is growing rather quickly with the 
growth of sC a t S C^10. 

Let us assess time to in the condition t^to of the applicability of Eqs. ( 1371) . ( 138|) and ( 1421 ). 
Let us rewrite to in accordance with Eqs. ( 1351) and ( |25l) as 



t ~R 2 /(dR 2 /dt). (43) 

Eq. ( 1431) shows that value represents the rate of relative increase of the area of the bubble 
surface at time t . 

For radius Ro introduced by means of Eq. ([6]), let us write 

i? ~40^, (44) 

where it is considered that the largest of the two values given in braces in Eq. ([6]) is, as a 
rule, value 2o~/IL In accordance with Eq. (l4~4"|) . relative deviation of pressure in the bubble 
of radius R^Ro from pressure II of the solution does not exceed 5%. Let us assume that 
II~10 6 ^r^l atm and cr~75 ^ (as in case of water). Then, from Eq. flU} we have 

#0-3 ■ 1(T 3 cm. (45) 

Let u s and the derivative in Eq. » by means of Fig. 2, using the valne that is 

typical for water: D~l,6 ■ 1CT 5 22- Further we shall obtain R 2 Q in Eq. (PD usine; the 
evaluation from Eq. (T45l) . The data on the dependence of the evaluation of time t Q on s( 
obtained in the described way are presented in Fig. 3. 
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FIG. 3: Dependence of the estimation for the time to ( s ) ° n the product s£. 



Using the analytical Eq. (I42p in Eq. (H3|) at the same value of diffusion coefficient D and 
the same evaluation from Eq. (H5|) . we have 

1,5- irr 1 



tc, 



s «>10). 



(46) 



Analytical Eq. (j4"B"l) and Fig. 3 provide exhaustive data on the evaluation of time t through- 
out all the area of s( > 0. 

Let us notice that the Eq. (136|) can be rewritten using Eq. (|35|) as 



R 2 (t)~R 2 t/t {m ), 



where values R and t evaluated previously are present. 
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(47) 



IV. STEADY THEORY CASE 



Let us study when 

«) 1/2 <1. (48) 

Thus, from Eq. (1521) it follows that 

a l ' 2 <l. (49) 

As we shall see further, inequality in Eq. (149]) will help us to solve a self-similar Eq. (129]) 
analytically. The same inequality will also stipulate the range of applicability of the steady 
theory. 

Let us begin by solving Eq. (129]) . Differentiating by b, one can make sure that the 
following asymptotic equation is observed: 

^ = 1 - (^)' 2 b 1/2 + ... (6 1 / 2 «1), (50) 

J i x \2J 

where the neglected components in terms of their value are even smaller than b. Using Eq. 
([50]) in a self-similar Eq. ( 1291 . considering only its main component, we shall have a = b, 
which, by virtue of Eq. (T49l) . confirms the condition 6 1//2 <d in Eq. (1501) . As a result, the 
solution of a self-similar Eq. (1291) . in case when Eq. (1491 is observed, in the principal order 
by a 1 / 2 looks as follows: 

b = a (a 1/2 <l). (51) 

Now let us define the range of applicability of the steady theory. Let us substitute Eq. 
fl5T]) in the self-similar Eq. (128]) . Thus we have 

n(p) =n QO + n g a^ ^ ex P ( _ ^r _ ^ + y) (a 1/2 <^l). (52) 

Let us prove that the concentration profile of the dissolved gas in Eq. ( |52l) . which was 
obtained at the limitation from above on a expressed in Eq. (T49|) . is steady in the area of 

Kp<2 1 / 2 /4a 1 / 2 (a 1 / 2 «l), (53) 

that exists simultaneously with limitation from Eq. (149]) . In the area stipulated by Eq. (153]) 
the absolute value in the exponent in Eq. (152]) is rather small (smaller than 1/16), and the 
exponent itself, consequently, is rather close to 1. Taking this fact into consideration and 
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also developing a in the multiplier in front of the integral in Eq. (|52l) . using definition from 
Eq. (130]) . in the area stipulated by Eq. fl53|) we shall with high precision deduce Eq. fl52|) to 



n(p)=n -^-^ (a^«l, l<p<2 1 / a /4a 1 / a ). (54) 

Eq. (|54|) corresponds to the steady concentration profile of the dissolved gas. We also 
provide limitations on a and on p; when those limitations are observed Eq. (154]) stems from 
a general self-similar Eq. (128]) for the nonsteady concentration profile of the dissolved gas. 

In order to find the rate of bubble growth, as it can be seen from Eq. (123]) . it is enough 
to know the derivative dn(p)/dp\ p= i, i.e. it is enough to know the concentration profile of 
the dissolved gas only in the infinitely narrow vicinity of the bubble. But even then one can 
obtain this data within the framework of the steady approximation of Eq. fl5^]) only when 
the limitation Eq. (H9]) is observed. Indeed, already when a 1//2 >2 1//2 /4 (when the limitation 
of Eq. (4.2) is no longer observed) the area stipulated by Eq. fl53]) does not exist at all; and, 
consequently, there is no such thing as an infinitely narrow vicinity of the bubble, where the 
concentration profile of the dissolved gas could be steady. 

We see that the limitation in Eq. ( fl9l is necessary for works aaaaaa based on 

the steady theory and devoted to the study of bubble growth. At the same time, it becomes 
clear that the rapid increase in the rate of bubble growth, which, as it was revealed in section 
ITT1 occurs at the increase of the product s(, and, consequently, at the increase of value a 
(when Eq. (4.2) is not observed) , is of a significantly nonsteady nature. 

Finally, let us apply the analytical solution from Eq. (15 ip of the self-similar Eq. ( 12 9 p to 
the problem of finding the rate of bubble growth. Substituting Eq. (15T]) in the self-similar 
Eq. (125]) . we have 

dR 2 /dt = 2Da (a 1/2 <l) (55) 

(for the sake of simplicity we do not point out the inequality t^to which is the condition for 
the observation of Eq. (125])). Exploiting Eq. ( 132]) in Eq. ( 155]) . we have 

dR 2 /dt = 2Ds( («) 1/2 <1). (56) 



Using the definition from Eq. ([3]), Eq. f]56]) can be rewritten as 

dR kTn^D 



dt RU 



C (57) 



Under the condition that Pa-CIT specified in the beginning of section [I] Eq. (157]) coincides 

n 

with Eq. (2.7) in [4j obtained at low gas solubility. 
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V. RANGE OF APPLICABILITY OF THE THEORY 

Let us first find the condition on the smallness of the deviation of the bubble temperature 
from the temperature of the solution. The transition of a molecule of gas from the solution 
to the bubble, as a rule, occurs with absorption of heat. As a result, the temperature of 
the bubble, which we will denote as T g , becomes smaller than temperature T of the solution 
(now we perceive T as the temperature of the solution at the infinite distance from the 
bubble) . In order to find the deviation T g — T, one preliminary has to find the nonsteady 
temperature profile T(r,t) around the bubble by solving the heat conductivity equation in 
the liquid solvent moving outside the bubble with the boundary condition T(r,t)\ r=OQ = T 
and at the intensity of absorption of heat by the bubble stipulated by the rate of its growth. 
This can be achieved, just as in section [Til by means of a self-similar method based on the 
substitution T(r,t) = T(p), where a self-similar variable p is given by the previously used 
equality f|T9|) . (Study of nonsteady heat conductivity problem is important for the growth 
of a liquid droplet in supersaturated vapor [9(). Using this method and obtaining T(p), we 
will have the following equality for the unknown deviation T g — T, where T g is defined as 



Here q is the heat of dissolution per one molecule of g rule, q > 0), k is the coefficient 

of heat conductivity of the liquid solvent in the diluted solution, 



X is the coefficient of thermal diffusivity of the liquid solvent in the diluted solution. For 
the rate of bubble growth in Eq. (1581) we have used Eq. fl25|) , which is enough to assert the 
smallness of heat effect, (but not enough to study it). 
Replacing value b with value eb in Eq. ( l40l) . we have 



T 9 = T(p)\ 



P =i- 




(58) 



e = D/x, 



(59) 




(60) 



From Eqs. (1581) and (1601 we have 





(61) 
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We shall assume values characteristic for water in ordinary conditions 8J: 

c cm 2 A erg , cm 2 

D~1,6-1(T 5 , k~6-10 4 ^— , x^l,4-10" 3 . 62 

s cm x s x K s 

Then, for value e introduced in Eq. ( 1591) we have 

£~1(T 2 . (63) 

As can be seen from Eqs. f )4T|) and fl63l) . limitation b>2 ■ 10 2 /e in Eq. flSTj) is equipotent 
to limitation a>10 2 which makes Eq. (|4T!) correct. Then, in accordance with Eqs. (ISTl) . 
flU and dHJ), we have 

T s _ r:s _±M% (a > 102 ). (64) 



Using Eq. (16^1) we have 
Let us accept that 



|T 9 -T| 10|g|L>n 9 



T kT 



,6 d y ne 



a (a>10 2 ). (65) 



n~10 b T„~300 K, \q\/kT~10. (66) 

cm 2 

We will use Eqs. fl62l) . fl66|) in Eq. fl65l) . and then expand n g by means of Eq. (TJJ. Thus we 
will have 

|T 9 -T|/T~l(T 4 a (a>10 2 ). (67) 

From Eq. fIBTl) it follows that 

(T.-TI/T^IO- 1 when a<10 3 , (68) 

where we do not refer to limitation a>10 2 in Eq. (1671) . as in accordance with physical meaning 
the value in the left part of Eq. (|68j) decreases with the decrease of value a. Inequality in 
Eq. fl68l) reveals the smallness of the relative deviation of the bubble temperature from the 
temperature of the solution even in the case of very large values of parameter a = s£<10 3 , 
where Eq.fl32l is observed. 

Now we will obtain the condition on the equalization of the temperature of the bubble. 
For time t g which is characteristic for this equalization we have the following estimate: 
tg^R 2 1 'x g , where Xg is the coefficient of thermal diffusivity of gas in the bubble. Let us now 
introduce characteristic time = R/(dR/dt) for the changing of the bubble radius. The 
fact that the temperature of the bubble manages to establish in the process of the bubble 
growth is equivalent to strong inequality t g <^itji, that is 

R/ Xg <£l/(dR/dt). (69) 
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By means of self-similar Eq. (124]) we shall rewrite Eq. (169]) as 



Db/x 9 ^l. (70) 

Using Eq. (HT1) . we have 

Z)6 6 .D o , 

= a 2 a>10 . 71 

From Eq. (1711) . taking Eq. (162]) and the fact that x 9 ~10 _1 — into consideration, it follows 
that 

Db/xg&O- 1 when a<l,8-10, (72) 

where we do not refer to limitation a>10 in Eq. (1711) . as value 6 decreases with the decrease 
of value a. Inequality in Eq. (1721) reveals the smallness of value _D6/x ff in the condition from 
Eq. ( 1701) at all values of sCi$l> 8 ' 10; where Eq.(l32l) is observed. 

Thus, if we consider conditions on the isothermal character of the bubble growth in Eqs. 
( 168]) and (175]) . the latter condition turns out to be stronger, that is 

<<1,8-10. (73) 

Still it leaves broad possibilities for the application of the isothermal theory. Indeed, the 
condition from Eq. (1731) in accordance with Eq. (142p admits even that 

< D~ 1 dR 2 /dt<l, 2 • 10 3 «<1, 8 • 10, (74) 

where the limitation s£>10 in Eq. (j4~2"]) is not mentioned as D~ 1 dR 2 /dt decreases with the 
decrease of s(. As for velocity 1/to of the relative increase of the area of the bubble surface 
at time to, the condition in Eq. ( 1731) in accordance with Eq. ( 1461) admits even that 

< l/t <2, 3 • 10 3 s- 1 «<1, 8-10), (75) 

where limitation s^lO in Eq. (146]) is not mentioned as l/to decreases with the decrease of 
s(. The ranges of alteration of values D~ l dR 2 /dt l/to (174"]) (1751) comprise three orders of 
magnitude. Therefore, the possibilities of the isothermal theory are rather broad. 

Let us prove that the observation of the condition from Eq. (175]) causes the observation 
of the condition on the diluteness of the solvent. In accordance with Eqs. (|4]), ([5]) and (ITBl) . 
we have 

n = ( s ( + s)U/kT. (76) 
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As s-ClO is correct in all cases, then from Eq. (I76j) in accordance with Eq. (I66I) it follows 
that 

n <4,3-10 20 cm- 3 «<1, 8 • 10). (77) 

The estimation from Eq. (177|) shows that, when the condition in Eq. (1731) is observed, the 
concentration no of the dissolved gas is approximately two orders less that the concentration 
n\ of the liquid solvent (n;~3 ■ 10 22 cm -3 for the case of water). This substantiates the 
diluteness of the solvent. 

Let us finally find the condition which is responsible for the observation of the mechanical 
equilibrium of the bubble and the solvent which was presupposed in the paper. We will 
use th to denote the time of mechanical relaxation of the bubble under the influence of 
the hydrodynamic interaction between the bubble and the solvent. As we consider large 
bubbles, i. e. cases when, in accordance with the estimation in Eq. ( 1451) . -R>3 • 10~ 3 cm, 
the following inequality is observed 

m Pp R 2 /Arj 2 > 1, (78) 

where p p and 7] are mass density and the viscosity of the solvent accordingly. Indeed, this 
fact is obvious from n~10 6 -^v, p„ = 1 and ri = 10 -2 g . As it has been shown 

Qcm z ' r P cm° 1 cmxs 

, when inequality in Eq. ( 1781) is observed, the mechanical relaxation of the bubble 

is of oscillation and attenuation character, while time th of this relaxation is given by the 

following equality: 

t h = Pp R 2 /2 v (79) 

n 

(Eq. (23) in [2]). Let us introduce characteristic time tjt = R/ (dR/dt) of the alteration of 
the radius of the bubble. The condition on the observation of the mechanical equilibrium of 
the bubble and the solvent is 

t h /tR<l- (80) 
Obtaining tn by means of Eq. (14"2"]1 and taking Eq, (1791) into account, we have 



From Eq. (JglJ), when£>~l, 6 • 10~ 5 ^f, p p = 1 and rj = 10 



<) 2 «>10). (81) 



t h /t R <5 -10- 1 «<1,8-10), (82) 
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where we do not point at limitation s£?^10 in Eq. (IHTj) . as in accordance with physical 
meaning value th/tR decreases with the decrease of value s(. In accordance with Eq. fl82l) . 
the inequality in Eq. (15U|) is observed in the vicinity of the boundary value even under the 
condition of Eq. f|75|) which, as it has been mentioned earlier provided extensive possibilities 
for the application of the theory. We will impose this condition on s( in Eq. (I82p in order 
to reveal it that the mechanical equilibrium between the bubble and the solvent and the 
equalization of temperature inside the bubble cease being observed at approximately the 
same values of s(. 

In conclusion we will note that the smallness of gas solubility was not exploited in the 
present paper. 
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